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ABSTRACT 

A new method is introduced allowing to solve exactly the reactions A+A— s>inert 
and A+A— s>A on the ID lattice with synchronous diffusional dynamics (simultaneous 
hopping of all particles). Exact connections are found relating densities and certain 
correlation properties of these two reactions at all times. Asymptotic behavior at large 
times as well as scaling form describing the regime of low initial density, are derived 
explicitly. 
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In this work we develop a new method of obtaining exact results for certain one- 
dimensional reaction-diffusion models. Dynamics of the ID Ising chain was first solved 
by Glauber [1]. This famous solution has lead to a host of results explored over the 
last three decades. However, recently there has been a resurgence of interest in the 
solvable ID models with both new emphases and new solution methods, developed 
by many authors [2-25]. The new emphases have been on those models which have 
no equilibrium states, such as irreversible reactions [2-17], cluster coarsening at phase 
separation [17-23], deposition processes [1,24], models of self-organized criticality [25], 
etc. The solution methods utilized a variety of different approaches both for continuous- 
time asynchronous, and to a lesser extent for synchronous (discrete time, simultaneous 
updating, cellular-automaton-type) dynamics [17,26]. 

While our approach can be applied to a larger class of models [27] , the present work 
is devoted to the diffusion-limited (i.e., instantaneous reaction on each encounter) 

coagulation, A-l-A— >A, and annihilation A-l-A— >inert, of particles synchronously 
hopping on the ID lattice. The models will be defined in detail later. Besides providing 
interesting examples of strongly nonclassical fluctuations, ID reactions also describe 
certain experimental systems [28] with the reactants being typically excitations in chain- 
like structures. 

Various results for these reactions have been reported in the literature [2-17] and 
they include several exact solutions, either in the continuum off-lattice limit [8-10,13] 
or on the lattice. Our results, besides providing previously unavailable (on-lattice) 
solutions for the synchronous case, resolve two long-standing theoretical issues in this 
field. Firstly, we elucidate and prove exactly the equivalence of the coagulation and 
annihilation reactions for all times and densities . This equivalence has been anticipated 
and explored by several authors [7-8,10,13-15] based on the asymptotic large-time results 
and on certain similarities between the correlation functions of both reactions. Here we 
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derive an explicit connection, relations (12)-(13) below. 

Secondly, earlier exact on-lattice calculations for A+A^inert, based on the map- 
ping [4] to the low-T dynamics of Ising interfaces [4,11-12,17], were for initial conditions 
corresponding to certain subtle correlations in the particle locations at time t = 0; see 
the discussion in [12]. The initially uncorrelated state was only possible for the (initial) 
density p = | (per site). As a result, the asymptotic large-time solution was not uni- 
versal. Our new results explicitly apply for the random initial distribution of arbitrary 
density p per site. The large-time behavior of both the coagulation and annihilation 
reactions is found to be universal and not dependent on p. The onset of this behavior 
occurs nonuniformly for small p. The latter regime can be described by a scaling form 
introduced in [29], to be reviewed later on, see relations (19) and (21) below. Our ex- 
act results confirm the proposed scaling and also yield explicit form of the appropriate 
scaling function. 

We consider a model of or 1 particles at lattice sites i of the ID lattice, at times 
t — 0, 1, 2, . . .. In the time step t ^ t + 1, all particles hop synchronously, to one of 
the neighboring sites i ± 1, with equal probability |. The decision which way to hop is 
done independently for each particle. However, if two particles end up at the same site 
they instantaneously "react" according to the annihilation rule (the particle number is 
reduced from 2 to 0) or coagulation rule (the number is reduced from 2 to 1). It is 
clear that this dynamics decouples the even-odd and odd-even lattices. Thus, we only 
consider particles at even lattice sites i = 0, ±2, ±4, ... at even times t = 0, 2, 4, . . ., and 
at odd lattice sites i = ±1, ±3, ... at odd times. At time t = all even lattice sites are 
randomly occupied (particle number 1) with probability p or empty (particle number 
0) with probability 1 — p. 

The method of exact solution employed here utilizes a system with a certain conser- 
vation property which allows simplification of the dynamical behavior. Thus, rather then 
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considering particles directly, we consider integer "charge" variables qi{t) = 0, 1,2, . . . 
at each lattice site (with the even-odd sublattice convention as before). The charges 
follow the dynamics introduced in the studies of self-organized criticality [25] : they ran- 
domly hop at each time step, similar to the particles. However, if two charges end up at 
the same site they stick together and move as a single charge (sum of the two original 
charges) at later time steps. Such dynamics can be represented by the rule 

q^{t + 1) = a,_i(t)Q,_i(t) + [1 - a,+i(t)] qi+i{t) , (1) 

where the hopping-decision variables aj{t) are or 1 with probability ^. Specifically, 
aj{t) — 1 indicates that the charge at j hopped to j -|- 1 in the time step t ^ t + 1, 
whereas the value corresponds to stepping to j — 1. 

The important feature of such dynamical rules [25] is that due to charge conserva- 
tion the sum of charges in r consecutive sites, 

Sr,j = Qj + qj+2 + • • • + gj+2(r-l) , (2) 

evolves with only two random decisions at the end-points needed, 

Sr,j{t+1) = aj-i{t)qj-i{t) + qj+i{t) + . . . + qj+2r-3{t) + [l - aj+2r-l{t)] qj+2r-l{t) , (3) 

where the random variables a at sites j — 1 and j + 2r — 1 enter. Relation (3) follows 
directly from (1) and (2). 

Consider now a function F{s) defined on the allowed s values 0,1,2, The 

averages 

fr{t)^{F{SrAm, (4) 

where the averaging is both over the random hopping decisions ai{t) and over trans- 
lationally invariant initial conditions, are translationally invariant, i.e., they do not 
depend on j. By using relation (3), one can easily verify that 

frit + 1) = ^ [/r+1 W + 2/.(t) + fr-l{t)] , (5) 
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where we define fo{t) = F{0) to have (5) apply for all r > 1. Earlier studies used 
this formulation with the choice F{s) = e*'^* and a further complication of allowing for 
added charge "fed in" at each time step, by including inhomogeneous term in (1) with 
values drawn from a fixed, time-independent distribution. The resulting characteristic 
functions fr{uj,t) can be used to study the self-similar steady state and dynamics [27] 
of charge systems with various charge variables, — real, positive, integer, etc., — in the 
framework of self-organized criticality. 

Our approach differs in that we employ functions F[s) which are not smooth. 
Specifically, let us define the empty interval indicator function. 



We relate the particle and charge systems as follows: a continuous interval of r sites 
is empty of particles and of charge simultaneously. Thus, we associate charge q > ^ 
at a given site with particle number 1, and charge g = with particle number 0. The 
initial state of or 1 particles at each site is represented by gi(0) = for sites i empty 
of particles, and by placing an arbitrary initial charge gi(0) > at sites occupied by 
particles. Then the dynamics of charges at the later times t > will also yield the correct 
dynamics of coagulating particles in the reaction A+A^A, which we will denote by the 
subscript A for brevity in what follows. 

The initial values of the averages are 



while the general difference 1 — frit) yields the probability that an interval of r consec- 
utive sites is empty. Empty interval probabilities have been considered in other studies 
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(6) 



/,(0) = (l-pr , 



(7) 



and for the particle density per site we get 



(8) 



of the coagulation reaction, e.g., [13]. Note that the system of recursion relations (5) for 
the coagulation reaction must be solved with the initial conditions (7) and "boundary 
condition" 



Our emphasis here will be on the density, c{t), which will be calculated exactly. 
However, let us first introduce the appropriate formulation for the annihilation reaction 
A+A— >inert, which will be denoted by the subscript for brevity. The idea to relate 
the aggregating "charge" system to 

hopping-particle models has been considered in [7,30]. Specifically, Spouge [7] used 
it to derive various exact continuum-limit expressions for both types of reaction (with 
asynchronous dynamics). Besides technical differences, the implementation here is more 
powerful than earlier variants. Firstly, synchronous dynamics and the associated charge- 
conservation property (3) allow derivation of explicit, discrete-time exact results in a 
rather straightforward manner. Secondly, extensions and generalizations are more easily 
identified [27]. 

For the annihilation reaction, we consider the even-occupancy indicator function. 



Indeed, the mapping to the particle system is now defined by associating even charges 
with empty sites and odd charges with sites occupied by particles. Initially, each empty 
site (particle number 0) is assigned an arbitrary even charge value q'(O), whereas every 
occupied site (particle number 1) is assigned an odd charge value. At later times t > 
the dynamics of charges will then describe the annihilation reaction A-l-A— >inert. 
Specifically, (t) = 1 — /i(t), while the difi'erence 1 — fr{t) denotes the probability that 
an r-interval has an even number of particles (0, 2, 4, ... , 2 [r/2] integer part) in it. 



fo{t) = F(0) = 1 . 



(9) 




(10) 
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The initial conditions for the annihilation case are more complicated than (7). 
Indeed, with a little bit of combinatorics one can check that the correct relation is 
fr{0) = [1 + (1 — 2py] jl. The boundary condition (9) is unchanged. This suggests 
consideration of the modified functions 

gr{t)^2fr{t)-l. (11) 

It is obvious that grif) satisfy, for with the initial density p^, the relations (5), (7) 
and (9) exactly identical to those satisfied by the original averages /r(t) for the reaction 
A provided we put 

"A = ^"0 ■ (12) 

Thus we discover that for the particular synchronous dynamics selected for the 
annihilation and coagulation reactions, the probabilities of finding even occupancy for 
the former are related to the probabilities of finding empty interval for the latter, with 
twice the initial density in coagulation as compared to annihilation. The precise relations 
can be quantified via the identity \grif)\Q = [/r(^)]^) with (12). For instance for the 
densities we find the relation 

CA(*) = 2c0(t), (13) 

which for this particular dynamics is exact for all times t = 1,2,3,... provided the 
random initial densities satisfy the relation (16) at t = 0; see (12). As mentioned earlier, 
several authors have explored [7-8,10,13-15] the asymptotic variants of this result for 
other dynamical rules. 

We now consider the coagulation reaction only and derive exact results for the den- 
sity by solving the recursions (5), with (7) and (9), by the generating function method. 
Only the outline of the actual solution steps will be presented. The emphasis will be 
on discussion of results for the density. Thus, we consider the generating functions 
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(I)r{u) = Y^t^o fr{t)u^ which Satisfy the relations (l)o{u) = 1/(1 — tt), 0r-(O) = (1 — pY, 
and 

u (0^+1 + 20^ + (Pr-i) = 4 [0^ - (1 - p)''] . (14) 

The solution of (14) is obtained by first eliminating the exponential-in-r dependence 
by defining (p^ = — pY'4'ri which yields an autonomous but inhomogeneous difference 
equation for 

V'r- However, the (constant) inhomogeneous term is then removed by the shift, 
ij}^ = 9r -\- Q- Finally, the resulting second-order difference equation for 9r is solved by 
the exponential form. Or = A'^Oq, where only one of the two roots of the characteristic 
equation gives the physically acceptable (i.e., regular at tt = 0) A value. The resulting 
expression is 

1 - (1 - w)i/2] u-^ [1/(1 -u) -&] + {!- pYQ , (15) 
e-' = l-u{2-pY/[A{l-p)]. (16) 

Let us now analyze the generating function for the density, 

l/(l-«)-0i(w) = ^c(t)«* = 2(l-«)-^/2 h+(2-p)p-^(l-«)^/2l"' _ (^7) 

t=o 

The final expression was obtained from (15) with (16) and required a rather cumbersome 
algebraic manipulation. 

Its form immediately illustrates that the large-time density, 

c{t ^ oo) = 2/V7rt , (18) 

is independent of p, as expected by the universality considerations referred to in the 
opening discussion; the residue of the leading singular term oc (1 — u)"-^/^ does not 
depend on p. 
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We also note that the onset of the large-time behavior in (17) is nonuniform in the 
limit p ^ 0. Indeed, it has been argued in the literature [29] that the universal large- 
time behavior sets in after the initial state is well "mixed" by the particle diffusion. The 
(dimensional) time scale for such a mixing is of order [(dimensional) density] "^/^/"D, 
where D is the dimensionality of space and V is the diffusion constant of the particle 
hopping. Here we have D= 1 and, in our dimensionless units, V =0(1). Thus, the 
onset of the large-time behavior in the limit of small initial densities will be described 
by a scaling form [29] which in our case is 

c{t ^ oo, p ^ 0) = t-^/^R{pH) . (19) 

The prefactor was selected to have the correct asymptotic form as t — > oo for fixed 
p > 0, with R{x = oo) = 2/V7r, while for short times, t <^ O (p~^), the reaction is 
ineffective: R{x — > 0) = y/x. 

In order to verify these expectations we need an explicit expression of the power 
series coefficients of (17). These can be obtained in terms of a hypergeometric function, 
or as polynomials in p, of degree t + 1, ior integer t values, or as a quadrature [27]. We 
favor the latter expression, 

4(l-p) 

_ p(2-p)^^+H2t + l)!! ^ z^dz 



(1 - p)*+l23*+2t! J (1 - ^)l/2 ' 



because it provides an explicit analytic continuation to all real t > values (even though 
the original model was only defined for integer t). 

Furthermore, the form (20) was already arranged to allow a relatively simple deriva- 
tion of the scaling relation (19). We find 

oo 
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which satisfies all the limiting properties as expected. Note that all the expressions for 
the coagulation reaction are meaningful for < p < 1. However, the relation to the 
annihilation reaction, (12), (13), etc., suggests that the formulas must be well defined 
for all < p < 2. Indeed, the appropriate results for the annihilation reaction are 
obtained by replacing p by 2p in all the "coagulation" expressions (and of course adding 
various other factors such as the coefficient \ in "translating" the density expression). 
Examination of (20) reveals that the powers of 1 — p in the denominator are cancelled 
by the appropriate factors from the integral (note 1 — p in the upper bound) so that the 
resulting expressions are indeed smooth and well defined near p = 1. 

In summary, we obtained new exact results, as well as the low-density scaling 
relation, for the ID coagulation and annihilation reactions with synchronous dynamics. 
These reactions are related to each other for all times. The large-time expressions 
are universal. We remind the reader that the notion of universality refers not only 
to the absence of dependence of results such as (18) on the initial density but also 
to the expectation that universal results apply to a class of models differing only in 
microscopic details of the dynamics, synchronous or asynchronous. However, for latter 
comparison the ID densities must be expressed per unit length, the dimensionless time 
related to the actual physical time r, etc. For instance the result (18) for the actual, 
dimensional density of the coagulation reaction then reduces to {2'kT>t)~^/'^ , etc. Our 
explicit results for the density of both reactions are universal when thus compared with 
other calculations available in the literature [2-17]. 

The author wishes to acknowledge the hospitality and financial assistance of the 
Institute for Theoretical Physics at the Technion, Israel Institute of Technology, where 
this work was completed. 
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